A simple model of universe with a polytropic equation of state 
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We construct a simple model of universe with a generalized equation of state p — (a + fcp 1//n )pc 2 
having a linear component p = ape 2 and a polytropic component p = £;p 1+1// "c 2 . For a = 1/3, 
n = 1 and fc = — 4/(3pp), where pp = 5.16 10 99 g/m 3 is the Planck density, this equation of state 
provides a model of the early universe without singularity describing the transition between the 
pre-radiation era and the radiation era. The universe starts from t = — oo but, when t < 0, its size 
is less than the Planck length lp = 1.62 10 -35 m. The universe undergoes an inflationary expansion 
that brings it to a size a± — 2.61 10 -6 m on a timescale of a few Planck times tp = 5.39 10 -44 s. 
When t S> tp, the universe decelerates and enters in the radiation era. For a = 0, n = — 1 
and k = — pA, where pA = 7.02 10 -24 g/m 3 is the cosmological density, this equation of state 
describes the transition from a decelerating universe dominated by baryonic and dark matter to an 
accelerating universe dominated by dark energy (second inflation). The transition takes place at a 
size a2 = 8.95 10 25 m corresponding to a time of the order of the cosmological time tA = 1.46 10 18 s. 
The present universe turns out to be just at the transition (to ~ tA)- This polytropic model reveals 
a nice "symmetry" between the early and late evolution of the universe, the cosmological constant A 
in the late universe playing a role similar to the Planck constant h in the early universe. We interpret 
the cosmological constant as a fundamental constant of nature describing the "cosmophysics" just 
like the Planck constant describes the microphysics. The Planck density and the cosmological 
density represent fundamental upper and lower bounds differing by 122 orders of magnitude. The 
cosmological constant "problem" may be a false problem. 



According to contemporary cosmology, the present en- 
ergy content of the universe is composed of approxi- 
mately 5% baryonic matter, 20% dark matter and 75% 
dark energy [l| . The expansion of the universe began in a 
tremendous inflationary burst driven by vacuum energy. 
Between 10 -35 and 10~ 33 seconds after the Big Bang, the 
universe expanded by a factor 10 30 @. The universe then 
entered in the radiation era and, when the temperature 
cooled down below approximately 10 3 K, in the matter 
era @. At present, it undergoes an accelerated expan- 
sion presumably due to the cosmological constant or to 
some form of dark energy Q . This corresponds to a sec- 
ond period of inflation. Despite the success of the stan- 
dard model, the nature of dark matter, dark energy, and 
the evolution of the very early universe (pre-radiation 
era), remains very mysterious and leads to many specu- 
lations. In this Letter, we propose to parameterize this 
indetcrmination by the concept of an equation of state. 
We show that the early and the late evolution of the 
universe can be described "symmetrically" by two poly- 
tropic equations of state with index n = 1 and n = — 1, 
respectively. Polytropic equations of state play an impor- 
tant role in astrophysics 0, § , statistical physics Q , and 
mathematical biology @ , and they may also be useful in 
cosmology. Our results are consistent with the standard 
ACDM model but refine it by removing the primordial 
singularity. 

We assume that the universe is isotropic and homoge- 
neous at large scales and contains a uniform perfect fluid 
of energy density e(t) = p(t)c 2 and pressure p(t). We also 
assume that the universe is flat in agreement with obser- 



vations of the cosmic microwave background (CMB) [l| 
It that case, the Einstein equations reduce to 
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where a(t) is the scale factor ("radius" of the universe) 
and H = a/ a is the Hubble parameter. These are the 
well-known Friedmann equations describing a non-static 
universe Q. The first equation can be viewed as an 
"equation of continuity" . For a given barotropic equa- 
tion of state p = p{p), it determines the relation between 
the density and the scale factor. Then, the temporal evo- 
lution of the scale factor is given by Eq. @. We have 
not written the cosmological constant A because its effect 
will be taken into account in the equation of state. Wc 
will also need the thermodynamical equation 



dp 1 , 2 . n 
dT = T {pC +P) > 



(3) 



which can be derived from the first principle of thermo- 
dynamics For a given barotropic equation of state 
p = p(p), this equation can be integrated to obtain the 
relation T = T(p) between the temperature and the den- 
sity. To close the system of equations, we have to specify 
an equation of state. We consider a generalized equation 
of state of the form 



p = (ap + hp 
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This is the sum of a standard linear equation of state 
p = ape 2 describing a universe dominated by radiation 
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(a = 1/3) or prcssureless matter (a = 0), and a poly- 
tropic equation of state p = kp 7 c 2 , where k is the poly- 
tropic constant and 7 = 1 + 1/n the polytropic index. An 
exhaustive study of this equation of state, considering all 
possible cases, is given in Q. In this Letter, we describe 
two specific models of particular physical interest. 

We first consider the case a = 1/3, n = 1 and k = 
— 4/(3/9p), which provides a model of the early universe 
without singularity. The equation of state can be 
rewritten as 



1 



V= oC 1 ~ 4 p/pp)p c 



(5) 



The linear part describes the radiation and the quadratic 
part may be due to Bose-Einstein condensates (BECs) 
with attractive self- interaction [loj . or have another ori- 
gin. The continuity equation ([T]) can be integrated into 



Pp 
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where a\ is a constant of integration. 

Pre-radiation era: When a -C a\, the density is ap- 
proximately constant: p ~ pp. Since this solution is 
expected to describe the very early universe, it is natu- 
ral to identify the constant pp with the Planck density 
which may represent a fundamental upper bound for the 
density. A constant value of the density gives rise to a 
phase of early inflation. From the Friedmann equation 
([2]), we find that the scale factor increases exponentially 
rapidly with time as 



a(t) ~ Ipe^W 2 *'** 



We have defined the time t = such that a(0) is equal 
to the Planck length lp. This universe exists at any time 
in the past (a — > and p — > pp for t — > —00), so there is 
no primordial singularity. 

Radiation era: When a 3> <Zi, we get pj pp ~ (a\ja) . 
When the density is "low" , the equation of state (|5|) re- 
duces to p = pc 2 /3 which describes the radiation era. 
The conservation of p ra d Q 4 implies that ppa\ — p ra d.o &o 
where p ra d,o is the present density of radiation and 
a = c/Hq = 1.32 10 26 m the present distance of cos- 
mological horizon determined by the Hubble constant 
H = 2.271CT 18 s~ 1 . Writing p ra dX) = &rad,opo where 
pa = 9.20 1 - 24 g/m 3 is the present density of the uni- 
verse given by Eq. ((2]), and £l r ad,o = 8.48 10~ 5 , we obtain 



ai = 2.61 10~ 6 m. 



(8) 



This scale is intermediate between the Planck length lp 
and the present size of the universe ao (ai/lp = 1.61 10 29 
and ai/a = 1.97 10~ 32 ). It gives the typical size of the 
universe at the end of the inflationary phase (or at the 
beginning of the radiation era). When a>oi, the Fried- 
mann equation ([2]) yields aja\ ~ (32n/3) 1 ^ 4 (t/tp) 1 ^ 2 and 
p/pp ~ (3/327r)(t P /t) 2 . 



The general solution: The general solution of the Fried- 
mann equation ([2j with Eq. ([6]) is 



s/{a/a x y + 1 - In 
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where C ~ —134 is a constant of integration determined 
such that a ~ lp at t = 0. The universe is accelerating 
when a < a c (i.e. p > p c ) and decelerating when a > a c 
(i.e. p < p c ) where p c = pp/2 and a c = 01. The time 
t c at which the universe starts decelerating is given by 
t c = (3/327r) 1 / 2 [ v / 2-ln(l + V2)-C]tp. This corresponds 
to the time at which the curve a(t) presents an inflexion 
point. It turns out that this inflexion point coincides with 
ai so it also marks the end of the inflation (t c = t\). 

Evolution of the temperature: The thermodynamical 
equation ([3]) with Eq. ([5]) can be integrated into 



T = Tp 
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T = T P 
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where Tp is a constant of integration. In the radiation 
era p < p P , Eq. dTUJ) reduces to p/p P - (tt 2 /15)(T/T p ) 4 
which is the Stefan-Boltzmann law provided that Tp is 
identified with the Planck temperature T P = 1.42 10 32 K. 
The temperature starts from T = at t = —00, in- 
creases exponentially rapidly during the inflation, reaches 
a maximum value T e , and decreases algebraically as 
T/Tp - (45/327T 3 ) 1 / 4 (i P /t) 1 / 2 during the radiation era. 
The point corresponding to the maximum temperature is 
Pe = pp/8, a e = 7^a u T e = (7/8) 7 / 4 (15/8Tr 2 ) 1 / 4 Tp. It 
is reached at a time t e = (3/32tt) 1 / 2 [\/8 - ln(l + + 
ln7/2 -C]t P . 

Early universe: The universe starts at t = — 00 with 
a vanishing radius a = 0, a finite density p = pp — 
5.16 10" g/m 3 , and a vanishing temperature T = 0. 
The universe exists at any time in the past and does 
not present any singularity. We define the "original" 
time t = ti = as the time at which the radius 
of the universe is equal to the Planck length. Thus 
ai = lp = 1.6210~ 35 m. The corresponding density 
and temperature are pi ~ pp = 5.16 10" g/m 3 and 
T t = 3.91 10- 205 T P = 5.5410^ 173 K. We note that 
quantum mechanics regularizes the finite time singular- 
ity present in the standard Big Bang theory. This is 
similar to finite size effects in second order phase tran- 
sitions (the Big Bang theory is recovered for fi = 0) 
0. We also note that the universe is very cold at 
t = ti = 0, unlike what is predicted by the Big Bang the- 
ory (a naive extrapolation of the law T oc t -1 / 2 leads to 
T(0) = +00). The universe first undergoes a phase of in- 
flation during which its radius and temperature increase 
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FIG. 1. Evolution of the scale factor a as a function of time. 
There is an early inflation (due to the Planck density pp) 
connecting the pre-radiation era to the radiation era. During 
the early inflation, the scale factor increases by 29 orders of 
magnitude in less than 10 -42 s. There is a late inflation (due 
to the cosmological density pa) connecting the matter era to 
the dark energy era. A present, we live just at the transition 
(bullet). 



FIG. 3. Evolution of the temperature T as a function of 
time. Before the early inflation, the universe is extremely 
cold (T < 1CT 173 K). During the early inflation, the tem- 
perature increases by 204 orders of magnitude in less than 
10 -42 s. During the radiation era, the temperature decreases 
algebraically. At present, T ~ 2.7 K (bullet). 



150 



100 



50 



o 

— 



-50 



1 1 1 

- s. 

— N.- 1 ? 
Pp \. 




Early inflation's^ 


Radiation era 




Matter era 

\ Pa 




~t A x » . 

a o Late inflation 

(dark energy) _ 


l l 


I 



-60 -50 -40 -30 -20 -10 10 20 30 40 50 60 
Log[t(s)] 

FIG. 2. Evolution of the density p as a function of time. It 
varies between two bounds pp and p\, fixed by fundamental 
constants of physics, producing two phases of inflation. 



exponentially rapidly while its density remains approx- 
imately constant. The inflation "starts" at t j = and 
ends at t\ = 23.3 tp = 1.251CP 42 s. During this very 
short lapse of time, the radius of the universe grows from 
a t = l P = 1.62 1(T 35 mtooi = 2.61 1CT 6 m, and the tem- 
perature grows from T z = 3.91 1CT 205 T P = 5.54 lCT 173 K 
to Ti = 0.278 T P = 3.93 10 31 K. By contrast, the den- 
sity does not change significatively: It goes from pi ~ 
p P = 5.1610"g/m 3 to Pl = 0.5 p P = 2.58 10" g/m 3 . 
After the inflation, the universe enters in the radiation 
era and, from that point, we recover the standard model 
Q. The radius increases algebraically as a oc t 1 / 2 while 
the density and the temperature decrease algebraically 
as p oc t~ 2 and T oc t^ 1 / 2 . The temperature achieves 
its maximum value T e = 0.523 T P = 7.40 10 31 K at 



t = t e = 23.6 t P = 1.2710~ 42 s. At that moment, the 
density is p e = 0.125 pp = 6.44 10 98 g/m 3 and the ra- 
dius a e = 1.63 ai = 4.2410~ 6 m. During the inflation, 
the universe is accelerating and during the radiation era 
it is decelerating. The transition takes place at a time 
t c = ti = 23.3 tp = 1.2510 _42 s coinciding with the end 
of the inflation (a c = ai). The evolution of the scale fac- 
tor, density and temperature as a function of time are 
represented in Figs. [T][3]in logarithmic scales. 

Scalar field: The phase of inflation in the very early 
universe is usually described by a scalar field @ . We have 
determined the ordinary scalar field (minimally coupled 
to gravity) leading to the equation of state (J5J) and found 
that it corresponds to a potential 



1 2 cosh 2 ?/. + 2 

V{ip) = -ppc -j— . 

o cosh tp 



(12) 



We now consider the case a = 0, n = — 1 and k = —p\, 
which provides a model of the late universe dominated 
by dark matter and dark energy. The equation of state 
can be rewritten as 



P 



-PAC 



(13) 



The pressure has a constant negative value. It belongs 
to the family of the generalized Chaplygin model p = 
-A/p a with a > -1 (here a = 0) [!,[ll|]. The continuity 
equation (Q~J) can be integrated into 



P = PA 



(14) 



where 0,2 is a constant of integration. 

Dark energy era: When a 3> a 2 , the density is ap- 
proximately constant: p ~ pa- Since this solution is 
expected to describe the very late universe, it is natural 
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to identify the constant pa with the cosmological den- 
sity p A = A/8ttG = 7.02 10~ 24 g/m 3 (dark energy) which 
may represent a fundamental lower bound for the density. 
A constant value of the density gives rise to a phase of late 
inflation. It is convenient to define a cosmological time 
*A = l/(Gp A ) 1/2 = (8tt/A) 1 /2 = 1.4610 18 s and a cos- 
mological length l A = ct A = (8ttc 2 /A) 1 / 2 = 4.38 10 26 m. 
These are the counterparts of the Planck scales for the 
late universe. From the Fricdmann equation ([2]), we find 
that the scale factor increases exponentially rapidly with 
time as 



a(t) ~i AC (8»/3) 1/a (*-*/)/*A, 



(15) 



with t f = 0.708t A = 1.0310 18 s (see below). This ex- 
ponentialgrowth corresponds to the so-called de Sitter 
solution [l(. 

Matter era: When a <C a 2 , we get p/p A ~ (a 2 /a) 3 . 
Because of the smallness of p A , the equation of state 
f| 1 3|) can be approximated by p ~ which describes the 



pressureless matter era @ . The conservation of p m a" 
implies that p A a 2 = p m ,o a o- Using p A = 51a,oPo an d 
p m .o = Q m ,oPo with f^A.o = 0.763 (dark energy) and 
&m,o = ^B.o + Qdai.o = 0.237 (baryons and dark mat- 
ter), we obtain 02/00 = (^m.o/^A.o) 1 ^ 3 hence 



a 2 = 8.9510 25 m. 



(16) 



This can be rewritten as 02 = 0.204/a- When a <C 02, the 
Friedmann equation yields a/a 2 ~ (67r) 1//3 (t/t A ) 2 / 3 
and p/p\ ~ (l/6Tr)(t A /t) 2 (Einstein-de Sitter solution). 

The general solution: The general solution of the Fried- 
mann equation ((2|) with Eq. (fl4|) is 



a 
a 2 



sinh 



2/3 



'67T- 



(17) 



Its asymptotic behavior gives Eq. (fT5j) with t/ = 
(3/87r) 1 / 2 [(2/3)ln2 + ln(l A /a 2 )]t A . The universe is de- 
celerating when a < a' c (i.e. p > p' c ) and accelerating 
when a > a' c (i.e. p < p' c ) where a' c = (l/2) 1 / 3 a 2 and 
p' c = 3p A . The time t' c at which the universe starts accel- 
erating is given by t' c = (l/67r) 1 / 2 argsinh(l/v / 2)tA- This 
corresponds to the time at which the curve a(t) presents 
a second inflexion point. 

Late universe: In the matter era, the radius increases 
algebraically as a oc t 2 ^ 3 while the density decreases 
algebraically as p oc t~ 2 (EdS). When a 3> a 2 , the 
universe enters in the dark energy era. It undergoes 
a late inflation (de Sitter) during which its radius in- 
creases exponentially rapidly while its density remains 
constant and equal to the cosmological density p A . The 
transition takes place at a 2 = 0.204^a = 8.95 1 25 m, 
92 = 2pa = 1.4010- 23 g/m 3 , t 2 = 0.203i A = 2.9710 17 s. 
In the matter era, the universe is decelerating and in the 
dark energy era it is accelerating. The time at which the 
universe starts accelerating is t c — 0.152£a = 2.22 10 17 s, 
corresponding to a radius a c = 0.162Za = 7.11 10 25 m and 
a density p c = ip A = 2.11 10" 23 g/m 3 . 



Present universe: The present size of the universe 
oq = 0.302/a = 1-32 1 26 m is precisely of the order of the 
scale a 2 = 0.204/ A = 8.95 1 25 m (o = 1.48o 2 ). There- 
fore, we live just at the transition between the matter 
era and the dark energy era (see bullets in Figs. Q] and 
[3]). The present density of the universe is po = 1.31pA = 
9.20 10 _24 g/m 3 and the present value of the deceleration 
parameter is qo = (1 — 30 A ,o)/2 = —0.645. The age of 
the universe is to = (l/67r) 1 / 2 argsinh[(ao/a2) 3 / 2 ]tA- Nu- 
merically, t = 0.310U = 4.54 10 17 s. 

Scalar field: In alternative theories to the cosmolog- 
ical constant, the present-day acceleration of the uni- 
verse (dark energy) is described by a scalar field called 
quintessence Q . We have determined the potential of the 
quintessence field associated with the equation of state 
(JTB"]) and found that it is given by 

V(f/>) = ip A c 2 (cosh 2 i> + 1), (V> > 0). (18) 

On the other hand, since w = p/ pc 2 = —p A /p is between 
— 1 and (wo — — ^a,o = —0.763), the equation of state 
(flB"]) can also be associated with a tachyon field Q. We 
have found that it corresponds to a potential 



COS i/j' 



(0 < ip < tt/2). 



(19) 



Conclusion: We have considered two simple polytropic 
equations of state with negative pressure and index n = 1 
and n = —1. These equations of state are relatively 
"symmetric" . They describe respectively a phase of infla- 
tion in the early universe (due to the Planck density pp) 
and in the late universe (due to the cosmological pa)- 
They also describe, in a unified manner, the transition 
from the pre-radiation era to the radiation era in the early 
universe, and the transition from the matter era to the 
dark energy era in the late universe. In the early universe, 
this amounts to summing the inverse Planck density and 
the inverse radiation density (1/p = 1/ pp + l/p ra d); in 
the late universe, this amounts to summing the matter 
density and the cosmological density (p = p m + pa). The 
whole evolution of the universe can be described by a 
single equation 



j + ^m,o [(a /a) 3 + {a /a 2 y 



Hq y (a/a ) 4 + (01/00) 

(20) 

combining the merits of these two models (the asterisk 
means that the matter term must be introduced for a 3> 
ai in order to avoid a spurious divergence when a — > 0). 
This equation can be written equivalently as 



H / Qrad.O 

Ho ~ V (a/a ) 4 + ( ai /a ) A (a/a ) 3 



+ O A , , (21) 



For ai = 0, we recover the standard ACDM model which 
presents a singularity at t = [l|. For a\ 7^ 0, the initial 
singularity is removed and the universe exists eternally 
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in past and future (aioniotic universe). This model is 
the most natural non-singular symmetric solution of the 
cosmological Einstein equations Q . 

Cosmological constant problem: It is oftentimes ar- 
gued that the dark energy density p\ represents the vac- 
uum energy density. Since the vacuum energy density 
is expected to be of the order of the Planck density 
Pp ~ 10 122 pa, this leads to the so-called cosmological 
constant problem [T2]. Actually, as illustrated in Fig. [2] 
the Planck density and the cosmological density repre- 
sent fundamental upper and lower density bounds acting 



in the early and late universe, respectively. Therefore, it 
is not surprising that they differ by 122 orders of mag- 
nitude. We view the cosmological constant as a funda- 
mental constant of nature describing the "cosmophysics" 
(late universe) in the same sense that the Planck constant 
describes the microphysics (early universe) . Accordingly, 
the origin of the dark energy density p\ should not be 
sought in quantum mechanics, but in pure general rela- 
tivity. In this sense, the cosmological constant "problem" 
may be a false problem. Of course, the origin of the cos- 
mological constant still remains to be understood. 
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